Quark number scaling of Vi in transverse kinetic energy and it's implications for 

coalescence models 
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We find that a simple extension of the coalescence model is sufficient to incorporate the perfect 
quark number scaling behavior of the elliptic flow in transverse kinetic energy, recently discovered by 
the PHENIX Collaboration. The flavor dependence of the elliptic flow can be consistently described 
in the low and intermediate pr if the transverse kinetic energy is conserved in the 2 — > 1 or 3 —* 1 
parton coalescence process at the hadronization. Thus suggesting the quark coalescence as a possible 
hadronization mechanism at low pr as well. 



PACS numbers: 25.75.Dw 

Quark recombination or coalesce models [l], H, 0, BL Q 
have beenproposed to explain the large baryon to meson 
ratio @, 0] and the number of constituent quark scal- 
ing (NCQ-scaling) of the hadron v 2 in pr 0,11] seen in 
RHIC data. Since coalescence models explicitly assume 
partonic degree of freedom, their apparent, albeit quali- 
tative, success in describing the data at intermediate px 
is a strong support of the idea that the partonic matter 
has been created in Au-Au collisions at RHIC. 

Coalescence models calculate the meson production 
via the convolution of Wiener functions for constituent 
quarks and the meson 0, El, 0, [H : 



distributions are driven by the space-time evolution of 
the partonic matter, which is typically described by hy- 
drodynamics or parton cascade. Coalescence models ex- 
plicitly assume that Fm can be factorized into the sin- 
gle quark distribution function, Fm — f(xi,pi)f(x2,p 2 ), 
with identical anisotropic phase space distribution (only 
V2 is considered here.) 



d 3 N 
dp 3 

U) 

where F is the joint q and q phase space distribution, 
Wm is the meson Wigner function, which is typically 
approximated by 

W M = $m{xi - x 2 ,Pi -P2)5{pt - Pt,i - PT.2) ( 2 ) 

The S function enforce the momentum conservation. It is 
straightforward to generalize these formula for baryons. 

There are several implicit assumptions in coalescence 
models. Only valence quark degrees of freedom are as- 
sumed in the coalescence calculation, and the dynamic 
gluon contribution (appears as higher Fock state) @ 
could lead to a 20% correction to the NCQ-scaling. There 
is a unitarity problem at low px @, H| due the quadratical 
dependence of the meson yield on the number of quarks, 
although this can be reconciled when one take into a fi- 
nite freeze out time and no quarks are used in more than 
one hadron (lfj |. In addition, coalescence process reduce 
the entropy. This can be partially resolved by including 
resonance decay of unstable hadrons [Til ]. 

As EqQ] indicates, the level of NCQ-scaling depends 
on both the final state parton phase space distribution, 
described by Fm, an d the hadronization process which 
is encoded in Wm- The constituent quark phase space 



f(x,p) = n(x,p) (1 + 2v 2 , q (x,p)cos(2cf))) 



(3) 



For a naive coalescence model where q and q are assumed 
to be co-moving : $m oc 5 3 (xi — x 2 )S 3 (pi — p 2 ), the 
spatially averaged v 2 for quark and meson is LU 



J d 3 Xi d 3 piF M {xi, pi, X2, P2 ) Wm (p, Pi ,P2, #1,^2) 
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where (A(x)) x = J d 3 xA(x). Because quark density 
n(x,p) and parton anisotropy v 2tq (x,p) could depend 
on quark local coordinate and momentum, NCQ-scaling 
in general does not follow as pointed in pjl, The 
space momentum correlation of the parton system before 
hadronization requires detailed study of the dynamical 
evolution of the system [l^, [H| • 

The level of NCQ-scaling also depends on the mod- 
elling of hadronization process. It is known that the fi- 
nite width of the hadron internal wave function can lead 
to 10-20% violation of the v 2 NCQ-scaling In the 

extreme case where wave function is very wide, NCQ- 
scaling can be destroyed completely [lq |. Since coales- 
cence deals with 2, 3 — > 1 processes, in general quarks or 
the coalescenced hadron need to be off mass shell [17j . 
and the interactions with the surrounding medium are 
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FIG. 1: (Color online) Figure taken from pl| a) A compila- 
tion of 1)2 /n as function of pr/n. b) The ratio to a common 
polynomial function. 



necessary to neutralize their virtuality. Thus the con- 
stituent quarks should be treated as quasi-particles and 
there is no reason a priori to assume the momentum con- 
servation in coalescence process (EqEJ. Our goal is to 
study the consequences of modifying this assumption. 

If coalescence is a relevant hadronization mechanism, 
it should be valid at all momentum. However, from the 
data point of view, NCQ-scaling of the «a, when plotted 
in p T , is not perfect (l8| . i.e. v 2 {pT/ l n)/n is not a univer- 
sal curve. This is illustrated by FigQ] where the ratios 
of the scaled elliptic flow (v 2 /n) for various hadrons to a 
combined polynomial fit are plotted as function of scaled 
transverse momentum (pr/n). NCQ-scaling only works 
at 20% level at pr > 2 GeV/c, then it breaks badly 
at pt < 2 GeV/c. The breaking has a systematic de- 
pendence on the mass. It undershoots the v 2 values for 
light mesons such as pion and overshoots the i> 2 for heavy 
baryons such as 5. Attempts have been made to reconcile 
the breaking of pion v 2 by including the resonance effect 
in the coalescence models [HI, [l9| . The agreement with 
data can be achieved qualitatively. We are not aware of 
any mechanisms within existing coalescence models that 
can accommodate the deviation of heavy baryons. 

Recently it was point out by PHENIX Collabora- 
tion [20I ] that the constituent quark scaling of i> 2 can be 
preserved at low pr if hadron transverse kinetic energy 
KEt = iut — m instead of px is used as a scaling vari- 
able. We duplicated the PHENIX finding in FigH The 
deviation from the universal v 2 (KET/n)/n scaling curve 
is very small, less than 10% in all cases. Motivated by 
this observation, we modify the momentum conservation 
relation into kinetic energy conservation relation in the 
coalescence formula Eq[2j 

Wm = $m(%i— %2,Pi— P2)S(rriT— mr,i— mr,2-m+mi+mj] 

(6) 
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FIG. 2: (Color online) Figure taken from |20( ] : a) A compila- 
tion of «2 /n as function of KEt /n. b) The ratio to a common 
polynomial function. 



Without very detailed calculations, we can make the fol- 
lowing remarks. 

1 . For the naive coalescence scenario, where only the 
co-moving q and q recombine into hadron, assuming 
spatially uniform quark density n(x,p) — n(p) or flow 
anisotropy V2, q (x 1 p) — V2,q(p), we obtain the usual NCQ- 
scaling relation for meson and baryon 0, @ , except it is 
in KEt, instead of px- 



v 2 ,m(KE t ) 
w 2 ,b(KE t ) = 
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3v 2 , q (— -) (8) 



We have confirmed these relations through a simple 
Monte-Carlo simulation of the naive coalescence process, 
where q or q are assumed to have mass of 350 MeV and an 
exponential spectra in mr- In addition, similar to what 
is done in other recombination models, the transverse 
momentum of the combining quarks in x or y direction 
are required to be within 0.24 GeV/c of each other. In 
essence the coalescence formula EqlT] becomes 



d 3 N 
dp 3 



d pid p 2 e 



-(mT,i+ m 



T - 2)/A S(KE T - KE T ,i - KE T , 2 ) 



6(A P - \p x .i -p X)2 |)6(A p - \p y>1 -p y ^\) 



(9) 



where A = 0.17 GeV and A p = 0.24 GeV/c is the cut off 
of relative momentum between the two quarks. Since we 
focus on mid-rapidity, longitudinal momentum direction 
has been ignored in this study. Results of our calculations 
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FIG. 3: (Color online) a) The input quark n 2 (KE T ). b) The 
ratio of the scaled V2 (u2(KEx/n)/n) for meson and baryon 
over quark i>2 (solid lines) and ratios of the scaled vi from 
Eq[7]and Eq[8]over quark V2 (dashed lines). 

are shown in Fig[3J The calculations are also compared 
with the analytical results given by Eq[7|and[8]in dashed 
lines. Clearly, the breaking of scaling due to the finite 
momentum spread in coalescence become significant only 
at KEt < 1 GeV. Above 1 GeV, both curves approach 
the value indicated by EqJT] and [8] (indicated by dashed 
lines). However, the split between meson and baryon in 
KEx/n would not be discernable given the precision of 
the data shown by FigH We noted that Eq[7] and Eq[S] 
are valid as long as the constituent quark spectra are 
exponential in mx, independent of the quark mass. 

2. The kinetic energy conservation condition is more 
than a simple improvement over the previous coales- 
cence models only at low px- It represents a quanti- 
tatively different recombination condition at all px- At 
large (pt 3> m), coalescence condition EqJ5] becomes 
Pt — Pt.i — Pt.2 = to — mi — 7712- It is different from mo- 
mentum conservation, because in general m ^ mi + 777,2. 
This can be clearly seen by comparing the bottom pan- 
els of FigfTjand Figf2] Going from px/n to KEx/n, one 
effectively shift high pt part of light meson V2 /n to the 
right and the low high px part of heavy baryon t;2 /n to 
the left by an amount on the order of the hadron mass. 
However, the resulting scaled vi curves falls on top of 
each other. 

3. In hydrodynamic model with Cooper- Fryer freeze 



out [2l|, hadrons are statistically produced in the local 
co-moving frame with certain radial flow. At low px, the 
elliptic flow should be proportional to p T according to 
the a general analytical expansion [22]. The Buda-Lund 
hydro model [H[ suggests a universal elliptic flow scaling 
: V2 = I\(w)/ Io{w),w oc p T /2rn,T) which, at lowpT (thus 
small w), reduce to v-i ~ 0.5w oc p T /2m = KEt- We also 
arrive this relation through a schematic derivation [241 ]. 

For any variable X and a linear relation between el- 
liptic flow and X : V2 oc X, NCQ-scaling in X is a 
trivial additive artifact: Vifn oc X/n. It was real- 
ized by [l3l [lij that either hydrodynamical mass order- 
ing or NCQ-scaling of V2 can results from quark coales- 
cence. In our point of view, the coalescence condition 
EqJS] provide a natural reconciliation between the sta- 
tistical hadronization and hadronization via coalescence. 
Since t>2 is proportional to KEt at low pt, «2 for various 
hadrons naturally fall on a universal curve in small KEt 
independent of whether they are scaled by the number of 
constituent quark or not. 

Due to the non-linear relation between pt and KEt, 
at low or intermediate px region where mass effect is 
important, NCQ-scaling can not hold in both pt and 
KEt simultaneously. At pt 3> m, a schematic hydro 
calculation gives v-i ~ tanh(£px) L 2l|, where £ depends on 
flow velocity. However, this relation gives neither NCQ- 
scaling in pr nor NCQ-scaling in KEt. 

A complete model should be able to describe not only 
elliptic flow but also the particle ratio, single particle 
spectra, etc. Our study is rather limited in this regard. 
We do not attempt to describe any observables other than 
elliptic flow. Various hadronic effects after quark coales- 
cence, such as the hadronic re-scattering and resonance 
decay, are ignored in current study. Our discussion is 
focused on the physical consequences when quark coales- 
cence is extended to low px- We have demonstrated that 
results from a simple coalescence approach agree with the 
data fairly well in KEt- The deviation of the pion v-i /n 
from pr/n scaling is naturally resolved when plotted in 
KE-r/n variable, without the need to rely on the reso- 
nance effect. The fact that RHIC data indicate strong 
flavor dependence of the V2 and strong violation of the 
NCQ-scaling of the v 2 in px suggests that KEt scaling of 
the ?J2 at intermediate px is a rather non-trivial observa- 
tion. It suggests that the quark coalescence might be an 
important hadronization mechanism at low px (pt < 2 
GeV/c). 

We thank M. Issah for pointing us to the Fig. 2 in his 
talk [23]. Discussions with R. Fries, C. M. Ko are appre- 
ciated. 
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For a thermal source and at pr >C tuvt, where vt is the 
radial flow velocity, 



j t (itlt — p~t ■ vt) ± 
f(p,x) ocexp( - ) 

7im T , , m(v(4>)) 2 
oc exp( — ) cc exp( \^ ) 



(10) 



2T 



we have ignored the <p dependence of the T and fi. For 
this simple scenario, the elliptic flow would be, 



V2 
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8mT 



(11) 



) y — {v) x \/{v) is the velocity anisotropy 
and should to first order be proportional to spacial 
anisotropy [2l]]: e v = Ae x + 0(pt). 
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